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It is shown that the boost-invariant and cylindrically non-symmetric hydrodynamic equations 
for baryon-free matter may be reduced to only two coupled differential equations. In the case 
where the system exhibits the cross-over phase transition, the standard numerical methods may be 
applied to solve these equations and the proposed scheme allows for a very convenient analysis of 
the cylindrically non-symmetric hydrodynamic expansion. 
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I. INTRODUCTION 



The success of the hadron production models based 
on the relativistic hydrodynamics 0, IS El IE IS 0] m 
describing the RHIC data brings the evidence that the 
hot and dense matterproduced at RHIC behaves like an 
almost perfect fluid yj. This situation, in turn, triggers 
interest in the further development of the hydrodynamic 
models, especially in the more detailed studies of the dis- 
sipative effects |E EE Ell- In this paper we consider 
the hydrodynamics of the perfect fluid and show that 
the boost-invariant and cylindrically non-symmetric rel- 
ativistic hydrodynamic equations for baryon-free matter 
may be comfortably reduced to only two coupled differen- 
tial equations. The effects of the cross-over phase transi- 
tion may be included in this scheme by the use of the tem- 
perature dependent sound velocity c s (T). As long as the 
function c s (T) satisfies the stability condition against the 
shock formation, the resulting equations may be solved 
with the help of standard numerical methods and the 
proposed scheme allows for a very convenient analysis 
of the cylindrically non-symmetric hydrodynamic expan- 
sion. The proposed method may be used to describe 
the evolution of matter produced in the central region of 
ultra-relativistic heavy-ion collisions, such as investigated 
in the present RHIC or future LHC experiments. 

The presented formalism is a direct generalization 
of the approach introduced by Baym, Friman, Blaizot, 
Soyeur, and Czyz [12]. In Ref. [T2] the boost-invariant 
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and cylindrically symmetric systems were considered, 
and the numerical calculations were performed only for 
the case of the constant sound velocity, c 2 s = ^. The ef- 
fects of the temperature dependent sound velocity were 
included in Ref. [13], where also the possibility of the 
non-zero initial transverse flow was considered in setting 
up the initial conditions for the hydrodynamic equations. 
In this paper we further relax the assumptions of Ref. 01 
by considering the cylindrically non-symmetric systems. 



II. HYDRODYNAMIC EQUATIONS 

The relativistic hydrodynamic equations of the perfect 
fluid follow from the energy-momentum conservation law 
and the assumption of the local equilibrium. For baryon- 
free matter they read 



vPd^Tu") = d v T, 
d^su") = 0, 



(1) 
(2) 



where T is the temperature, s is the entropy density, and 
= 7(1, v) is the four-velocity of the fluid. Due to 
the normalization of the four-velocity, only three out of 
four equations appearing in JJJ are independent. The as- 
sumption of the boost-invariance introduces another con- 
straint, hence, in this case Eqs. JTj and @ are reduced 
to three independent equations [lj| 
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FIG. 1: Decomposition of the flow velocity vector in the plane 
z — 0. In our approach we use the magnitude of the flow v 
and the angle a as two independent quantities. 



Here t, r = yj x 1 + y 2 , and (j) = t&n~ 1 (y/x) are the 
time and space coordinates which parameterize the plane 
z = (due to the boost-invariance of the system, the val- 
ues of all physical quantities at z ^ may be obtained 

directly by the Lorentz boosts). The quantity v is the 

_ i 

magnitude of the fluid velocity, 7 = (l — v 2 ) 2 is the 
corresponding Lorentz factor, while a is the function de- 
scribing direction of the flow, a = tan -1 (vt / 'v r) , see 
Fig. [I] The differential operator ^ in © is defined by 
the formula 



FIG. 2: Sound velocity squared shown as a function of the 
temperature (with the critical temperature T c — 170 MeV). 
The presented form is obtained by matchi ng t he lattice results 
with the resonance gas model calculation The cross- 

over phase transition is considered, hence no sudden decrease 
(dip) of the sound velocity is observed at T — T c . 



The condition against the shock formation has the form 
E03 



1 2 , rpdCa d (sc 8 \ 
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We note that this condition is naturally fulfilled in our 
case, since the derivative of c s (T) is always positive. 
With the help of the substitutions: 
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Equations © are three equations for four unknown 
functions: T, s, v, and a. To form a closed system of 
equations they should be supplemented by the equation 
of state, i.e., by the relation connecting T and s. Alter- 
natively, the equation of state may be included by the 
use of the temperature dependent sound velocity 



){T) 



dP _ s dT 
~de ~ T~ds' 



(5) 



and by the use of the potential <£(T) defined by the equa- 
tion 



(6) 



The form of the function c 2 (T) used in our calculations 
is shown in Fig. This form, with T c = 170 MeV, follows 
from matching the lattice results with the hadron reso- 
nance gas calculation pjj LL3] • We note that integration 
of Eq. (|JJ allows us to express in terms of the temper- 
ature, $ = $t(T), or to express temperature in terms of 
T = T$($). The functions $ T ( T ) and T$($), corre- 
sponding to the sound velocity function shown in Fig. EH 
are presented and discussed in more detail in Ref. 



tanh#, a± = exp($ ± 0), 



(8) 



the sum and the difference of the first two equations in 
© may be rewritten as 
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while the third equation in © is 
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Equations © and (fTU)) are three equations for three un- 
known functions: a+(t, r, (/>), a_(t, r, 0), and a(t, r, (j)). 
We note, that the velocity v and the potential are 
functions of a + and a . 



$=iln(a + a_). (11) 



Similarly, the sound velocity appearing in Eqs. © and 
(fTU|) may be also represented as a function of a + and a_, 



c s (T) 



ln(a + a_) 



(12) 



III. BOUNDARY CONDITIONS 

A. With cylindrical symmetry 

In the special case of the cylindrical symmetry all terms 
in Eq. (|T0)) vanish (9<£/90 = and a = 0, since we do 
not consider the possibility that the system rotates as a 
whole with a = tt/ 2), whereas Eqs. © are reduced to 
Eqs. (2.24) of Ref. O 



da± (v ±c s ) da± c s 
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a± = 0. (13) 



Due to the symmetry reasons, the velocity v should van- 
ish at r = 0, 



v(t,r = 0) = 0. 



(14) 



This condition is achieved by demanding that the two 
functions a+(t, r) and a_(£, r) may be expressed in terms 
of a single function a(t, r) with the help of the prescrip- 
tion 



a+(t, r) = a(t,r), r > 0, 
a_(£, r) = a(£, — r), r < 0. 



(15) 



The ansatz (|T5)> and the structure of Eqs. (fT3|) indicate 
that the equation for a_ (t, r) may be obtained from the 
equation for a+(t, r) by the substitution: r — > — r. This 
observation further means that two equations (|13|) may 
be reduced to a single equation for the function a(t, r) 
with the range of the variable r extended to negative 
values, 







da (v + c 3 ) da c s 
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(16) 



Furthermore, the ansatz (|T5)> and Eq. (|TT|) automatically 
yield the desired boundary condition for the temperature 



9T(t,r) 



dr 



(17) 



The scheme presented above is the basis of the approach 
used in Ref. [lj]. Below we develop the analogous scheme 
for cylindrically non-symmetric systems. 



B. Without cylindrical symmetry 

In this paper we consider the collisions of identical nu- 
clei with atomic number A which collide moving initially 
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FIG. 3: Construction of the functions a+(t, r, 0) and 
a~(t,r,<p) with the help of the function a(£,r, 0). The func- 
tion a(£, r, cj)) is symmetrically extended to negative values of 
r. 



along the z-axis. The impact parameter b points in the 
x-direction. The center of the first nucleus in the trans- 
verse plane is placed at xi = (xi, yi) = (—b/2, 0), and of 
the second nucleus at X2 = (#2,2/2) = (b/2, 0). In such a 
case, similarly to the symmetric case considered above, 
we also require that the magnitude of the flow vanishes 
at x = (x, y) = (0, 0), i.e., for r = we have 



v(t,r = 0,0) 0. 



(18) 



This condition is fulfilled naturally by the ansatz analo- 
gous to Eq. (p"5)) . namely 



a + (t,r,</>) = a(t,r,</>), r > 0, 
a_(£, r, 4>) = a(t, — r, 0), r < 0. 



(19) 



We supplement the ansatz (|T§)) by the definition of the 
function a{t, r, 0) for the negative arguments, see Fig. 03 



a(t, — r, 0) = a(t, r, < 



r > 0. 



(20) 



With the help of the definitions (fTTH) and pi. Eqs. © 
may be reduced to a single equation for the function 

a(t,r, 0), 
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a = 0, 



(21) 



Here, similarly to the cylindrically symmetric case, the 
range of the variable r is extended to negative values. 
Eq. (|2D should be solved together with Eq. (|TU)L where 
the range of r may be also extended trivially to negative 
values (definitions (fT§|) and (Q imply that this equa- 
tion is in fact invariant under transformation: r — > — r). 
The use of the polar coordinates in the transverse plane 
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requires also that all functions at (j) = and (j) = 2tt 
are equal: a(t, r, 0) = a(t, r, 2ir), a(t : r, 0) = a(t,r, 27r). 
We also note that Eqs. (fT§|) and (EQt yield the following 
boundary conditions for the temperature and the func- 
tion a(t, r, (j)), 



dT{t, r, ( 



dr 



0, 



<9a(t, r, ( 



7—0 



<9r 



0. 



(22) 



r=0 



In addition, multiplying the equation for a + in (0 by a_ 
and subtracting the equation for a_ multiplied by a + one 
finds 



1 dT(t, r, 



0. 



(23) 



r=0 



IV. INITIAL CONDITIONS 
A. Temperature 

In the following we assume that the hydrodynamic evo- 
lution starts at a typical time t = t = 1 fm. We assume 
also that the initial temperature profile is connected with 
the number of participating nucleons 

1/3 



T(to, x, y) = const 



dN p 
dx dy 



(24) 



where we use the formula Q 



dxdy \ 2 




In Eq. (Ji n = 40 mb is the inelastic nucleon-nucleon 
cross section and Ta (x, y) is the nucleus thickness func- 
tion 



T A (x,y) = / dzp(x,y,z) . 



(26) 



Here p{r) is the nuclear density profile given by the 
Woods-Saxon function with a conventional choice of the 
parameters used for the gold nucleus: 

0.17fm- 3 , 



a = 0.54 fm, A = 197. 



(27) 

The idea to use Eq. (EH follows from the assumption 
that the initially produced entropy density <j(to,x,y) is 
proportional to the number of the nucleons participating 
in the collision at the position (x,y). Since the consid- 
ered systems are initially very hot (with the temperature 
exceeding the critical temperature T c ), they may be con- 
sidered as systems of massless particles, where the en- 
tropy density is proportional to the third power of the 
temperature. In this way we arrive at Eq. (|21)k 



B. Velocity 



The initial form of the functions v(t, r, <p) and a(t, r, < 



is 



H r 



v(t ,r,<l>) = v (r)= - 7 ^- 

a(to,r,0) = 0, (28) 

where Ho is a parameter defining the initial transverse 
flow formed in the pre-equilibrium stage. In the present 
calculations we use a very small value Ho =0.001 fm -1 . 
We note that the finite value of the initial transverse flow 
makes all the terms on the right-hand-side of Eq. (|TU)) 
finite for r > 0. At r = these terms are also finite due 
to the boundary conditions (fT£|) and 

The two initial conditions, Eqs. (EH) and (Q ; may 
be included in the initial form of the function a(r) if we 
define 



/I + v 0( r ) 

a{t = t , r, 0) = a T (r, <j>) \ = n , (29) 



where 

ar(r, (j)) — exp < $t 



const 



dN p 
dx dy 



1/3 



(30) 



V. RESULTS 

The examples of our results are shown in Figs. 01 - 
El The three sets of plots correspond to three different 
values of the impact parameter; b = 3.9, 7.1, and 9.2 fm. 
These values are typical for the three centrality classes: 
c = - 20%, c = 20 - 40%, and c = 40 - 60%. The 
relation between the (average) impact parameter and the 
centrality class is obtained from the formula, 



1 



dcb(c) 



4r cUL 



3/2 



3 C n 



~5 (31) 



where the geometric relation between the centrality and 
the impact parameter is used, c = 6 2 /(4rg) [l3. In all 
considered cases the initial central temperature equals 
twice the critical temperature, To = T(to,0,0) = 2T C . 

The parts a) of Figs. 01- El show the temperature pro- 
files for different values of time: t = 1, 4, 7, 10, 13, 16 fm. 
The solid lines represent the temperature profiles along 
the x-axis (0 = 0), while the dashed lines represent the 
profiles along the y-axis (0 = 7r/2). One can observe that 
during the whole considered evolution time, 1 fm < t < 
16 fm, the ^/-extension of the system remains larger than 
the ^-extension. However, the relative magnitude of this 
effect decreases with time, indicating that the cylindri- 
cal symmetry is gradually restored as the time increases. 
In these parts of the plots one can also notice the effect 
of the phase transition; initially the system cools down 
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rather rapidly, later the cooling down is delayed and the 
main visible effect is the increase of the volume of the sys- 
tem. We note that this behavior is related to the sudden 
decrease of the sound velocity in the region T w T c . 

The parts b) of Figs. El-Elshow the isotherms in the t — 
r space, again for (j) = (solid lines) and (j) = tt/2 (dashed 
lines) . The pairs of isotherms indicate different values of 
the temperature. They start at T = 1.8 T c and go down 
to T = 0.2 T c , with a step of 0.2 T c . It is interesting to 
observe (especially for b = 7.1 fm and 9.2 fm) that the 
solid and dashed lines cross each other. This effect means 
that the central (relatively hotter) part of the system 
acquires a pumpkin-like shape during the evolution of 
the system. Such pumpkin-like regions, however, shrink 
and disappear during further expansion. 

In the parts c) of Figs. Ej-Elthe profiles of the func- 
tion a(t, r, (j>) are shown for t = 1, 4, 7, 10, 13, 16 fm and 
(j) = 7r/4. For = and (j) — ir/2 the function a(t, r, </>) 
vanishes due to the symmetry reasons. In the first and 
third quadrant the values of a are predominantly nega- 
tive, while in the second and fourth quadrant the values 
of a are positive. This behavior characterizes the direc- 
tion of the flow which has the tendency to change the ini- 
tial almond shape into a cylindrically symmetric shape. 
We have checked that positive values of a at (j) — 7r/4, 
seen at the edge of the system, are caused by the over- 
lapping tails of the Woods-Saxon distributions. 

In the parts d) the velocity profiles are shown, again 
for t = 1, 4, 10, 16 fm. Similarly to the previously dis- 
cussed figures, the solid lines are the profiles for = 
(x-direction) , whereas the dashed lines are the profiles for 
(j) = 7r/2 (^/-direction). One can observe that the magni- 
tude in the x-direction is larger than the magnitude in the 
^/-direction, which is an expected hydrodynamic behav- 
ior caused by larger pressure gradients in the x-direction. 
Exactly this effect is responsible for the observed az- 
imuthal asymmetry of the transverse-momentum spectra, 
quantified by the values of the V2 coefficient. 

Finally, in the parts e) we show the contour lines of 



the temperature, again for different values of time. These 
plots visualize the time development of the system. The 
arrows describe the magnitude and the direction of the 
flow (for better recognition the angle a is magnified by a 
factor of 3). 



VI. SUMMARY 

In this paper we introduced a new and concise 
treatment of the boost-invariant and cylindrically non- 
symmetric relativistic hydrodynamic equations. This ap- 
proach may be used to describe the evolution of matter 
produced in the central region of ultra-relativistic heavy- 
ion collisions, such as investigated in the present RHIC 
or future LHC experiments. The presented formalism 
is a direct generalization of the approach introduced by 
Baym et al. in Ref. [lj]. In the studied case, the symme- 
try of the problem allows us to rewrite the hydrodynamic 
equations as only two coupled differential equations, (fTU|) 
and (|2T|L which automatically lead to the fulfillment of 
the requested boundary conditions for the velocity and 
the temperature at the center of the system. The ef- 
fects of the phase transition are included in this scheme 
by the use of the temperature dependent sound veloc- 
ity. As long as the sound velocity satisfies the condition 
against the shock formation, the resulting hydrodynamic 
equations may be solved numerically with the help of the 
standard methods. We note that besides the sound ve- 
locity, no other thermodynamic variables are necessary 
for solving the hydrodynamic equations in this case. The 
presented results of the hydrodynamic calculations, sup- 
plemented with the appropriate freeze-out model, may 
be used to calculate different physical observables. In 
particular, the coefficient of the elliptic-flow, V2, may be 
extracted. To achieve this task, in the closest future we 
intend to combine our hydrodynamic approach with the 
statistical Monte-Carlo model THERMINAT0R [11]. 
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FIG. 4: Time development of matter characterized by the initial conditions ((^S|> - (p?T7|> with Ho = O.OOlfm 
T = T(t o ,0,0) = 2 T c , and b = 3.9 fm. 




FIG. 5: Same as Fig. Hbut with b = 7.1 fm. 




FIG. 6: Same as Fig. Hbut with b = 9.2 fm. 



